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Abstract 

We make use of the phase space density approach to discuss gravitino as a warm 
dark matter candidate. Barring fine tuning between the reheat temperature in the 
Universe and superparticle masses, we find that warm gravitinos have both appropriate 
total mass density, = ^dm — 0.2, and suitable primordial phase space density at 
low momenta provided that their mass is in the range 1 keV < rriQ < 15 keV, the 
reheat temperature in the Universe is low, Tr < 10 TeV, and masses of some of the 
superparticles are sufficiently small, M < 350 GeV. The latter property implies that 
the gravitino warm dark matter scenario will be either ruled out or supported by the 
LHC experiments. 

1 Introduction and summary 

The predictions of the ACDM model are in outstanding consistency with the bulk of cos- 
mological observations [1] (see also Ref. [2j and references therein). Yet there are clouds 
above the coUisionless cold dark matter scenario, which have to do with cosmic structure 
at subgalactic scales. Three most notable of them are missing satellites [3], cuspy galactic 
density profiles [1] and too low angular momenta of spiral galaxies [5J. All these suggest that 
CDM may be too cold, i.e. that the vanishing primordial velocity dispersion of dark matter 
particles may be problematic. Hence, one is naturally lead to consider warm dark matter 
(WDM) scenarios [HI [3 [H [9] . 

There are several ways to describe the difference between WDM and CDM scenarios. The 
most prominent one is that warm particles filter primordial power spectrum on small scales, 
and thus the formation of small halos is suppressed. The filtering scale must be small enough, 
since the power spectrum shows no significant deviations from the CDM prediction on scales 
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Figure 1: Linear matter power spectrum for standard ACDM cosmology (daslied line) and 
AWDM (solid lines) assuming the distribution of WDM particles as given by ([T]) with masses 
m = 1, 5, 10, 15, 20 and 30 keV and = 5mssm- 



within reach of current observations. This leads to constraints on the primordial velocity 
dispersion of WDM particles |ilOj. On the other hand, in order to improve on structure 
formation, the filtering scale must be of the order of the scale of missing satellites, which is 
believed to be of order 10^ - lO^M^ fTTl. 





We have calculated linear matter power spectrum in AWDM cosmology assuming that 
dark matter particles have the Fermi-Dirac primordial distribution function, normalized to 
correct present total density: 



67rC(3) mTo^g^^ eP/^o.-// + 1 ' ^ ' 

where m is the WDM particle mass, Tg^e// = ( ^ ) ] 9* and g^.^ = If are the effective 
number of relativistic degrees of freedom at the epoch of dark matter particle production and 
at present epoch, respectively. To this end we have modified the Boltzmann evolution equa- 
tions implemented in the Code for Anisotropics in the Microwave Background (camb) [12j. 
Figure [1] presents the resulting AWDM power spectrum for m = 1,5, 10, 15,20 and 30 keV 
(solid) in comparison with ACDM (dashed). In the WDM case, is chosen to be equal 



2 



to gussM = 228.75, the maximum number of relativistic degrees of freedom equilibrated in 
plasma in the framework of MSSM. One concludes that the power spectrum is suppressed by 
about an order of magnitude on the scales corresponding to IO^Mq and smaller provided the 
WDM particle mass is smaller than about 15 keV. Thus we consider a particle as a WDM 
candidate if its mass obeys 

m < 15 keV . (2) 

Of course, this is an indicative figure, not a strict upper limit. 

Another way to quantify the notion of warm dark matter is to make use of the phase 
space density approach [13], [HI [15] . Its key ingredient is the ratio between the mass density 
and the cube of the one-dimensional velocity dispersion in a given volume, Q = p/cr^. On 
the one hand, this quantity is measurable in galactic halos; on the other hand, it can be 
used as an estimator for coarse-grained distribution function of halo particles. Namely, for 
non-relativistic dark matter particles 

n 



Q m 



4 



(ip2)3/2 



where m is the mass of these particles and n is their average number density in a halo. 
Assuming that the coarse-grained distribution of halo particles is isotropic, fhaioip,''^) = 
fhaio{p,r), one estimates 
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[J fhaioip, r)d^p] 
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fhalo{P*,r) , 



where is a typical momentum of the dark matter particles. In this way the magnitude of 
the coarse-grained distribution function in galactic halos is estimated as 

f''''^" - 33/2 ^4 • ^^'> 

Coarse-grained distribution function is known to decrease during violent relaxation in 
collisionless systems [I6] . Hence, the primordial phase space density of dark matter particles 
cannot be lower than that observed in dark halos. This leads to the Tremaine-Gunn-like con- 
straints on dark matter models fl3[ [T5] . The strongest among these constraints are obtained 
by making use of the highest phase space densities observed in dark halos, namely those 
of dwarf spheroidal galaxies (dSph) [HI [H] . dSph's are the most dark matter dominated 
compact objects, and seem to be hosted by the smallest halos containing dark matter [11]. 
In recently discovered objects Coma Berenices, Leo IV and Canes Venaciti II, the value of 
Q ranges from 5 • 10~^ to 2 ■ 10~^ [T71 . In what follows we use the first, more 

^ ^ (km/s)"' (km/s)"* ' ' ' 

conservative value, 

Q = 5.10-^^. (4) 
(km/s) 
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By requiring that the primordial distribution function exceeds the coarse-grained one, / > 
fhaio, one arrives at the constraint 

33/2^V >Q- (5) 

This constraint gives rise to a reasonably well defined lower bound on m in a given model. 

If the primordial distribution is such that ([5]) is barely satisfied, the formation of high-Q 
objects like dSph's is suppressed. In fact, it may be suppressed even for larger /, since the 
coarse-grained distribution function may decrease considerably during the evolution. The 
parameter 

A = 

Q 

shows how strongly the coarse-grained distribution function / must be diluted due to relax- 
ation processes in order that the formation of dense compact dark matter halos be suppressed. 
It is known from simulations that the phase space density decreases during the structure for- 
mation. In particular, during the nonlinear stage it decreases by a factor of 10^ to 10^ [18j, 
or possibly higher. Hence, the primordial distribution function of WDM particles should be 
such that A > 10^ — 10^. At least naively, obtaining the dilution factor in a given model 
in the ballpark A = 1 — 10'^ would indicate that the primordial phase space density is just 
right to make dwarf galaxies but not even more compact objects. Interestingly, we will find 
that A is indeed in this ballpark for WDM gravitinos obeying ([2]). 

As discussed in Ref. [19], only a fraction of dark matter particles should definitely have 
high phase space density obeying (jS]). This fraction u is estimated as the relative contribution 
of dSph's into the total mass density of dark matter. Using the dwarf number density 
Udwarf — 7 ■ 10~^ Mpc~^ from Ref. [20] and assuming the average dSph mass of order 
lO^M© [11], one estimates 

u ^ ^ 10-5 . (6) 

"DM 

To be on conservative side, we impose the constraint ([5]) on this fraction of WDM particles 
only. Also, we calculate the value of A for this fraction. One expects that once the right 
fraction of the dark matter particles has the high phase space density, the most compact 
objects are produced in right numbers (and not overproduced). We have checked that our 
results would change very little if we used an estimate for i/ differing from even by an 
order of magnitude, i.e. u = 10^^ — 10^^. 

To summarize, we consider a dark matter model viable if a fraction ([6]) of its particles 
has primordial phase space density obeying ([5]) with Q given by Eq. (jll). 

In this paper we make use of this phase space density criterion together with the bound 
([2]) to examine light gravitino as a warm dark matter candidate, assuming that R-parity is 
conserved and hence gravitino is stable. We find that gravitino mass should be in the range 

lkeV<mG<15keV, 
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cf. i. 

In the early Universe, light gravitinos are produced in decays of superparticles and in 
scattering processes [211 1221 [231 121] • For so light gravitinos, their production in decays of 
superparticles plays an important role [22]. We consider this mechanism in Sec. 12.11 where 
we also evaluate the spectrum of produced gravitinos. In Sec. 12.21 we discuss gravitino 
production in scattering processes. The latter mechanism operates most efficiently at the 
highest possible temperatures in the early Universe, so the requirement that gravitinos are 
not overproduced restricts severely the reheat temperature Tr, cf. [251 126]; we find that Tr 
must be at most in the TeV range. 

Most notably, gravitinos serve as warm dark matter candidates only if other superparticles 
are rather light. We find that superparticles whose mass M is below the reheat temperature 
should obey 

M < 350 GeV , (7) 

otherwise gravitinos are overproduced in their decays and in scattering and/or relic gravitinos 
are too cold. Barring fine tuning between the reheat temperature in the Universe and 
superparticle masses, this means that gravitino as warm dark matter candidate will soon be 
either ruled out or supported by the LHC experiments. 

The bound ([7D is to be compared to the experimental bounds on masses of gluino and 
quarks of the 1st and 2nd generations, Mg ^ > 250 — 325 GeV [2]. Given the narrow interval 
between these bounds, we find it disfavored that squarks and gluinos participate in gravitino 
production processes. Hence, we elaborate also on a scenario with relatively light colorless 
superparticles whose masses M obey ([7]), heavy squarks and gluinos, and reheat temperature 
in between, 

M < « M^^g . (8) 

In this scenario, squarks and gluinos do not play any role in gravitino production, while 
the important production processes are decays and collisions of sleptons, charginos and 
neutralinos. We find that in this case, the overall picture is consistent in rather wide range 
of parameters, with the reheat temperature extending up to 10 TeV. 

It is worth noting that for light gravitino we consider in this paper, the lifetime of next- 
to-lightest superparticles (NLSP) is short, t^lsp ^ 2 ■ 10~^ s. Thus, their decays after 
decoupling are not hazardous for BBN. On the other hand, in the mass range of gravitino 
and superpartners we have found favored, one has tj^lsp ^ 5 ■ 10"'' s and thus the NLSP 
decay length (neglecting 7-factor) is in the range 

160 m < cr < 7 km . 

So, with gravitino WDM, it is likely that NLSP (if chargeless and colorless) will freely travel 
through the LHC detectors. 
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2 Gravitino production mechanisms 



2.1 Production in decays 

We begin with the study of the hght gravitino production in two-body decays of thermahzed 
superparticles, assuming that the reheat temperature in the Universe exceeds considerably 
the masses of these superparticles. Let us first find the distribution function of gravitinos 
produced in decays of one kind of superparticles with mass M. At time t, these superparticles 
have thermal distribution function fth{p,t). 

The Boltzmann equation for the gravitino distribution function f{p,t) is 

where H{t) is the Hubble parameter. In this section we consider the contribution into the 
collision term / that comes from two-body decays and is given by 

^ ^ "^'^ :{27rr6^'\P-p-p')UiP,t)\M\\ 



2|p| J 2E(27r)3 2|p'|(27r; 



3 ' 



Here P, p, p' are the 3-momenta of the decaying superparticle, gravitino and another decay 
product (SM particle), respectively; E = V + is the energy of the decaying particle. 
The amplitude A4 is related to the decay rate in the rest frame of the decaying particle as 
I Alp = IQttMT. Neglecting the masses of particles in the final state, one ha£| m 



We note that both decays and scattering processes produce longitudinal gravitino (gold- 
stino), so the number of gravitino helicity states effectively equals two. We also note that 
in the parameter range of interest, the inverse 2^1 processes, leading to disappearance of 
gravitino, have negligible rates. 

Upon integrating over the momentum of the SM particle and over the direction of P, the 
collision term takes the following form. 



Mr 

p2 



fth{P,t) dE 



^This formula, generally speaking, does not work in the Higgs-higgsino sector, and in some region of 
the parameter space the corresponding decays are suppressed. We treat higgsinos on equal footing with 
other charginos and neutralinos in what follows. Refining this approximation would not change our results 
considerably. 
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where 

_ M2 

Efnin P H Z 

is the minimum energy of the decaying particle capable of producing gravitino of momentum 
p. Of particular interest for what follows is the low momentum region, p <^ M, T. In that 

case E,nin ^ M, i.e., slow gravitinos arc born in peculiar decays of fast moving superparticlcs, 
which produce gravitinos in a narrow backward cone. For this reason, the efficient production 
of slow gravitinos occurs at temperatures T > M^/p ^ M. As we will see shortly, for 
relatively low reheat temperatures Tr this results in a non-trivial shape of the gravitino 
spectrum at low momenta, with a cutoff at p/T ~ M^/T^. 

It is convenient to take comoving momentum q = a{t) p as the argument of the gravitino 
distribution function. Here a{t) is the scale factor, whose present value is normalized to 
unity, a{to) — 1. The Boltzmann equation takes the form 



df{q,t) _ MT 
dt g2 

It can be easily integrated, giving 



oo 

a\t) J fthiP,t)dE. 



t oo 

/MV C 
dt'^a\t') j UP,t')dE, 



tR 

where tu refers to the beginning of the thermal phase of the cosmological evolution after 
reheating. Hereafter we assume that the production of gravitinos is negligible at the reheating 
epoch. This is of course an arbitrary assumption reflecting our ignorance of the reheating 
mechanism; we expect that the gravitino production at reheating, if any, would make the 
regions of favored gravitino and superparticle masses even narrower as compared to the 
regions presented below. 

Since the thermal distribution function fth{P, t) of the decaying particles depends on the 
ratio E/T{t) only, it is convenient to trade the integration over production time for the 
integration over temperature. To this end we use the entropy conservation and the relation 



T — y with Mpi = -^Piy 8^;^5 valid at the radiation domination epoch. Thus, if the 
gravitino distribution function had not been distorted by structure formation, at the present 
epoch it would have been given by 

y MVM*T^.f 7 fE\ 
mto)=J dT^^^ J dE. 
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Here Tn 



1/3 



; (7* and g^o = H are the effective number of relativistic degrees 



^0,e// = To 

of freedom at gravitino production and at present epoch, respectively; in the framework of 
MSSM with all superparticles relativistic in the plasma = (^mssm = 228.75 and Tq ^/j ~ 
0.7 K. 

Changing the variables {T,E) — > (^z = = and performing the integration over x 
we obtain finally the following result for the primordial distribution function expressed in 
terms of the momenta redshifted to the present epoch, 

2 



f{p)^f{q,to 



where 
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Tf 
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Tr 



37r3/2y^m|Mpi 
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M 



2Tj 



(10) 



with 



P 



+ 



Oe// 



Toe// 4J^ p 

The corresponding spectrum ^ = A-Kp^fip) for Tr ^ M is shown in the left panel of Fig. [2] 
in comparison with the thermal spectrum at temperature To,e// and the same total number 
of particles. It is seen that gravitinos produced in decays have lower average momentum 
and in this sense are cooler than thermal ones. The overall shape of the spectrum is not of 
particular interest for our purposes, however: the formation of compact objects like dSph's 
depends on the low-momentum part of the spectrum, where the phase space density f{p) is 
high. 

As we alluded to above, the low-momentum part of the spectrum depends in a peculiar 
way on the ratio of the mass of the decaying particle M to the reheat temperature Tr. This 
comes out from Eq. (fTOj) . in particular, through the lower limit of integration Zmin- For 
Tr — i> 00, the distribution function at low momenta, p -C To^eff, is given by 



fip) 



2\/5 fT^ 



37r3/2 
v^C(5/2) 



'0,e// 



p 



1/2 



z^'l^fthiz) dz 



gCde 



T 



0,eff 
P 



1/2 



where g is the number of helicity states of the decaying particle, and Cdec = 1 for bosons and 
Cdec 273) ^'^^ fermions. As expected, the combination m^f entering h 
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Figure 2: Left: Spectra of gravitinos produced in decays of thermalized fermions (solid) and 
bosons (dashed) in comparison with the Fermi-Dirac (dash-dotted) spectrum at temperature To^e//; 
normahzed to the same total number of particles. Right: Low momentum part of the distribution 
functions of gravitinos with itIq = 10 keV produced in the decays of bosons with M = 200 GeV 
for Tr oo (solid line) and Tr = 4M (dashed line). Dash-dotted line: Fermi-Dirac distribution 
at temperature To^eff, normalized to the same total number of particles. 

with the gravitino mass, m^f oc m^, so that only hght gravitinos are warm. We note in 
passing that the distribution function derived in this way is unbounded as p ^ 0. However, 
the correct distribution function of gravitinos — fermions with effectively two helicity states 
— cannot exceed the value 2/(27r)^ because of Pauli-blocking. To take this into account we 
simply cut the distribution function at 2/(27r)'^ wherever the calculated distribution function 
exceeds this value. In fact, this procedure is used almost nowhere in the parameter space we 
consider in this paper, as the calculated distribution function almost never exceeds 2/(27r)^. 

For finite Tr, but still Tr > M, the distribution function no longer peaks at p ^ 0. 
Instead, it has a rather broad peak at p/T ~ M'^/T^ and exponentially decays towards p 0. 
This is shown in the right panel of Fig. [2], where we also compare the distribution function of 
gravitinos produced in decays with the thermal distribution function at temperature To,e// 
normalized to the same total number of produced particles. We again see that the gravitinos 
produced in decays are substantially cooler than fermions with thermal distribution, as the 
maximum phase space density is substantially higher in the former case. 

It is clear from ([9]) that the largest contribution into the gravitino production comes from 
the heaviest superparticles that have ever been relativistic in cosmic plasma. To get an idea of 
numerics, let us consider the case in which gi, bosonic and Qf fermionic superparticle degrees 
of freedom have one and the same mass M, and the reheat temperature is substantially 
higher than M. Then the present number density of gravitinos produced in decays of these 
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superparticles is 



and the present mass density of these gravitinos is given by 

= !^ « 8 . 10- L + %,) f liEX) f^^X . (u) 

° Pc V 16 / V 9. y V "G / \ WOGeVj ^ > 

A crude estimate for the gravitino mass is obtained by assuming that the distribution function 
of = 10"^ of gravitinos is roughly comparable to the Pauli-blocking value, / = 2/(27r)^. 
Then the condition ([5]) corresponds to 

> 1 keV . (12) 

As an example, if the heaviest superparticles are squarks of the 1st and 2nd generations and 
gluinos, as motivated by mSUGRA, if they have the same mass and the reheat temperature 
is high enough so that these particles were relativistic in the cosmic plasma, then gb = Oq = 
4 • 3 ■ 4 = 48, gj = gg = 2 ■ 8 = 16 and g^ = gussM- Making use of the estimate (ITTl) 
and the upper limit on warm gravitino mass ([2]), we find in this example that the common 
mass of squarks and gluinos must be rather small, Mg^g < 350 GeV, otherwise gravitinos are 
overproduced. We will refine these estimates in Sec. [31 

2.2 Gravitino production in scattering 

Gravitino production in scattering processes has been worked out in Refs. [221 123] using the 
Braaten-Yuan prescription and hard loop resummation. It has been reconsidered recently in 
Ref. pi] with the results substantially different from those of Refs. [221 [23] in some regions of 
parameter space. We will use the approach of Refs. [221 [23] with understanding that there is 
considerable uncertainty both in gravitino production rate and in their spectrum, especially 
at relatively low temperatures, T ~ M. We will further comment on this uncertainty in 
Section [3l 

The contribution of scattering into the gravitino production is dominated by the processes 
involving the heaviest superparticles which have ever been relativistic in the cosmic plasma. 
Furthermore, this contribution strongly depends on whether or not these superparticles are 
colored. In what follows we consider two scenarios which we think are representative for 
realistic supersymmetric extensions of the Standard Model. Our analysis below is straight- 
forwardly redone for the general case, but given the unknown superparticle spectrum and 
the uncertainty in ([S]), considering these simple scenarios will be sufficient for our purposes. 
The first scenario has been described in the end of Sec. 12.11 in this scenario the heaviest are 
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squarks of the 1st and 2nd generations and gluinos, and we assume that they all have the 
same mass M and that the reheat temperature exceeds M . Given the experimental bounds, 
M > 250 — 325 GeV, it is clear already from the preliminary discussion in the end of Sec. 12.11 
that this scenario may be consistent only in a rather narrow range of the parameter space. 
Hence, we discuss also the second scenario, which is defined by the relation ([H]) where M is 
the common mass of sleptons, charginos and neutralinos. In the second scenario, squarks 
and gluinos play no role in the gravitino production in the early Universe. Given the strong 
dependence on the mass M, varying the rest of SUSY parameters in either scenario does not 
lead to significant changes of our results. 

For the first, squark-gluino scenario, the results of Ref. ^23j apply directly, so the mass 
density of gravitinos produced in scattering is given by 

ng «c..,.M„ C-i] f^^V (l^) (^] . (13) 

° \9sJ VlOOGeV; V '"G / VlTeV; ' ' 

where gs is the strong coupling constant at the energy scale Tr, and Ug = 0.732, kg = 1.271. 

For the second, color-singlet scenario, the results of Ref. [23] have to be modified. To this 
end, we consider electroweak scattering processes only and omit the contributions of reactions 
with external squarks. Also, we omit the squark contributions into the thermal masses of 
the gauge bosons. The overall gravitino production cross section depends on thermal masses 
rrith as In [T/mti,) and thus grows as the thermal mass decreases. As a result, the gravitino 
production cross section in our scenario is nearly 80% of the electroweak part obtained in 
Ref. [23], although 1/3 of all processes are omitted. Using the modified cross sections and 
g^: = 142.75, we find for the present mass density of gravitinos produced in 2 2 processes 
in primordial plasma: 



with modified constant factors Ua ~ (0.152, 0.372) and scales in logarithms ka ~ (1.52, 1.52). 
Here a = 1 and a = 2 refer to the gauge groups f/(l)y and SU{2)l, respectively, with the 
gauge couphngs ga = {g\g)- 

The estimates ([T51) and f[T^ have considerable uncertainties related to infrared problems 
existing in field theory at finite temperature. These will translate into uncertainties in our 
estimates presented in Section [31 

3 Results 

There are three criteria the WDM model with light gravitino should satisfy. First, as dis- 
cussed in Introduction, gravitino would serve as warm dark matter provided its mass satisfies 
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the upper bound ([2]), ttiq < 15 keV. Another criterion is that the present gravitino mass 
density should be equal to the observed dark matter density. In both scenarios of Sec. 12. 2[ 
the total gravitino mass density is the sum of contributions due to the decay and scattering 
processes, so that one requires + Q% = ^dm ~ 0.2. This requirement gives one relation 
between the three parameters, the masses itlq, M and reheat temperature Tr in each sce- 
nario. For the first, squark-gluino scenario we make use of f|T3l) as well as ffTTl) with = 48, 
Qf = IQ and = (^mssm- For the second, color-singlet scenario the appropriate expressions 
are ([H]) and ^ with = = 3 ■ (4 + 2) = 18, ^/ = = 4 ■ 2 + 2 ■ 4 = 16 and = 142.75. 
In either case, scanning the reheat temperature from M upwards, we observe from 

Eqs. flTTl) . f[T^ and f[T^ that this criterion gives a lower bound on the gravitino mass for 
given M. 

The third criterion is discussed in Sec. 1: about 10^^ of gravitinos should have the 
primordial phase space density obeying ([5]) with m = m^,. This criterion gives a lower 
bound on the gravitino mass for given M and Tr. This bound has to do with the magnitude 
of the gravitino distribution function at low momenta where this function is large. Instead 
of calculating the low momentum part of the distribution function of gravitinos produced in 
the scattering processes, we first use the lower bound on the overall distribution function, 
which is obtained by neglecting altogether the contribution of the scattering processes into 
the distribution function in the low momentum region. The lower bounds on itIq obtained 
within this decay dominance approximation are overestimated in comparison with those one 
would obtain by the complete treatment. To get an idea of the uncertainty introduced by 
approximating the low momentum part of the distribution function by the contribution of the 
decay processes only, we then add the contribution from scattering assuming that gravitinos 
produced in the latter way have thermal-shaped distribution ([T]), but normalized to the total 
mass density, Eqs. f|T3l) and f|T^ in the first and second scenario, respectively. 

Within either approximation, for each set of parameters we find the value / of the phase 
space density, such that 10~^ of gravitinos have the distribution function exceeding /, and 
require that / obeys ([5]) at an allowed point in the parameter space. 

The resulting bounds in (M, m^) plane are shown in Fig. [3] and in Fig. H] for the first and 
second scenario, respectively. Contours of equal Tr/M are plotted with dashed-dotted lines. 
At the same time, these contours correspond to constant fractions of gravitinos produced in 
scattering and decay channels, providing together the correct present mass density of dark 
matter, = Qdm- The labels show the fraction of gravitino produced in decays. The 
shaded regions are allowed by both ([2]) and ([5]). We also show the lines of equal dilution 
factor A; the dashed lines correspond to the decay dominance approximation, while the 
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Figure 3: Allowed region of masses (shaded) in a scenario with heavy gluinos and quarks of the 1st 
and 2nd generations, Mq = Mg = M, and Tr > M. Contours of equal dilution factor A are also 
shown (solid and dashed lines). The dashed lines correspond to the decay dominance approximation, 
while the solid lines are obtained under the assumption that the scattering contribution to the 
low momentum part of the distribution function has the form (jlSp . Contours of equal Tr/M 
are shown with dash-dotted lines, on which the fraction of gravitinos produced in decays is also 
indicated. Conservative experimental lower bound on masses of gluinos and squarks of the 1st and 
2nd generations is indicated by solid vertical line. 

solid lines are obtained under the assumption that the scattering contribution to the low 
momentum part of the distribution function has the form (fT5ll . 

In view of substantial uncertainty in the production of gravitinos in scattering, the esti- 
mates for the reheat temperature should be considered as indicative only. This is particularly 
relevant for the upper left parts of Figs. [3] and HI where production in scattering dominates 
over production in decays. Also, the estimates for the dilution factor A are uncertain in 
these parts of the parameter space, due to the large uncertainty in the low momentum part 
of the spectrum of gravitinos produced in scattering. This is reflected by the fact that dashed 
and solid lines deviate significantly from each other in the upper left parts of Figs. [3] and HI 
Furthermore, even though the thermal-shaped distribution (ITSl) is a plausible approximation. 
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Figure 4: Same as in Fig. [3l but for a scenario with color singlet superparticles of equal mass M, 
heavy squarks and gluinos, and intermediate reheat temperature, M < Tr <^ ^^,g- 



we cannot exclude the possibility that scattering contribution to the distribution function of 
gravitinos is much larger at low momenta as compared to f|T5|l . In the latter case the lines 
of equal dilution factor A would shift even further down. In any case, the most conservative 
lower bound on the gravitino mass independent of the distribution function is given by fll2p . 

Irrespectively of these uncertainties, we see that in both scenarios, the relevant superpar- 
ticle masses must be rather low, M < 320 — 350 GeV, provided that the reheat temperature 
is Tr > M. Extending the mass range of superparticles towards larger M in either scenario 
would require increasingly strong fine tuning between the reheat temperature and these 
masses. This fine tuning is needed to ensure that superparticles are non-relativistic and 
hence not so numerous, but have just right abundance at the beginning of the thermal stage 
of the cosmological evolution to produce just right number of gravitinos. We consider this 
possibility implausible. 

Figure [5] shows the same bounds as in Fig. H] in (Tr, m^) plane. On dash-dotted lines the 
total density of gravitinos produced in both channels is equal to the observed dark matter 
density for indicated superpartner masses M. 
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Figure 5: Region with gravitino WDM in the (TjijItIq) plane for the same scenario as in Fig. [H 
Contours of equal M are shown with dash-dotted lines. 

We conclude that unlike in the WIMP case, gravitino WDM does not automatically have 
the present mass density in the right ballpark. If the heaviest superparticles are squarks and 
gluinos, and they were relativistic in the cosmic plasma (the first scenario), the allowed range 
of parameters is rather narrow, as seen from Fig. [31 We consider least contrived the possibility 
that the masses of sleptons, charginos and neutralinos are in the range M = 150 — 300 GeV, 
the reheat temperature is Tr = 200 GeV — 10 TeV and the masses of gluinos and squarks 
are higher, Mg^g ^ Tr (second scenario). Then for masses rriQ = 1 — 15 keV, gravitinos can 
indeed serve as warm dark matter particles. In any case, gravitino as warm dark matter 
candidate will be either ruled out or supported by the LHC experiments. 
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